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We propose the quantitative mean-field theory of mechanical unfolding of a globule formed by 
long flexible homopolymer chain collapsed in poor solvent and subjected to extensional deformation. 
We demonstrate that depending on the degree of polymerization and solvent quality (quantified 
by the Flory-Huggins x parameter) the mechanical unfolding of the collapsed chain may either 
occur continuously (by passing a sequence of uniformly elongated configurations) or involves intra- 
molecular micro-phase coexistence of a collapsed and a stretched segment followed by an abrupt 
unravelling transition. The force-extension curves are obtained and quantitatively compared to 
our recent results of numerical self-consistent field (SCF) simulations. The phase diagrams for 
extended homopolymer chains in poor solvent comprising one- and two-phase regions are calculated 
for different chain length or/and solvent quality. 

I. INTRODUCTION 

Recent developments in polymer micromanipulation techniques such as atomic force microscopy (AFM) or optical 
tweezers have made it possible to subject an individual polymer molecule to mechanical deformation and to investigate 
the response of the molecule to this stimulus l^"^. The output provided by such single molecule force spectroscopy 
experiments is commonly expressed in the form of the force-extension curves. The ultimate goal of theory of mi- 
cromechanical deformation of macromolecules is to rationalize the shape of the force-extension curves obtained in 
experiments, i.e. to correlate different patterns (regimes) on the force-extension curves with specific conformational 
transitions that occur upon mechanical action on the macromolecule. 

The most complex patterns (non-monotonic, oscillating dependences) in the force-extension curves are observed 
upon mechanical manipulations with biological polymers, i.e., globular proteins, DNA and polysaccharides; these 
patterns are attributed to cascades of intra-molecular conformational transitions related to unfolding of hierarchically 
organized intra-molecular structures which are stabilized by strong intra-molecular attractive forces operating in 
aqueous media. However, even such a simple system as a homopolymer globule collapsed in poor solvent may 
exhibit a non-trivial behavior under extensional deformation. Analysis of this model system is highly important for 
rationalization of the relationships of mechanical unfolding of more complex globular structures and supra-molecular 
assemblies of biomacromolecules. 

There exist two modes of micromechanical manipulation depending on whether the value of deformation (extension) 
D or that of the force / applied to the macromolecule is chosen as the control parametei'^. From the point of view 
of the statistical mechanics, one can consider these two situations as corresponding to different statistical ensembles: 
The former case will be referred to as the constant extension ensemble, or Z?-ensemble, whereas the latter - as the 
constant force ensemble, or /-ensemble. 

The first theory of globule unfolding in the constant extension ensemble was developed by Halperin and ZhulingPl 
in the framework of the scaling approach. This theory is strictly applicable close to the theta-point since the solvent 
strength is accounted for in scaling terms and only for long chains comprising large number of thermal blobs. It was 
demonstrated that there exist three regimes of globule deformation: (i) At weak extensional deformations the globule 
has an elongated shape, Figure [l]b, the reaction force grows linearly with extension, (ii) The further increase in the 
end-to-end distance leads to appearance of a new intra-molecular micro-phase, i.e., a (strongly) stretched chain part. 
In a wide range of extensions, a microphase segregation takes place and the globule acquires a "tadpole" conformation 
with a prolate globular "head" and a stretched "tail" , Figure [l] b. The deformation is accompanied by progressive 
unfolding of the globular core and occurs at almost constant reaction force. This behavior has an analogy with the 
Rayleigh instability in a liquid dropletP under additional constraint of connectivity of solvophobic monomer units in 
the chain, (iii) At strong deformations the tadpole conformation becomes unstable and the globule completely unfolds 
to a uniformly extended chain. Figure [TJc. The reaction force grows again in this regime. According to the theory of 
Halperin and Zhulina, transitions between the regimes occur continuously upon an increase in the end-to-end distance 
D] in particular, the tadpole regime ends up when the size of the tadpole head becomes comparable to the thermal 
blob size (equal the thickness of the tadpole tail) or, in other words, when the head becomes indistinguishable from 
the tail. 

Later, Cooke and William^ considered stretching of a dry (solvent free) globule. They found that the tadpole 
conformation becomes unstable when the head of the tadpole still contains a large number of monomers, 7i ^ N^^^ 
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(in contrast to n ~ 1 that follows from prediction of the theory of Halperin and Zhulinaf^, if it is extrapolated to the 
dry globule case). This "unraveling transition"^ (i.e. complete globule unfolding) is accompanied by a sharp drop in 
the reaction force. 

Craig and Terentje\ii2l extended the Cooke and Williams approach to the case of collapsed semiflexible polymers and 
took into account chain stiffness and finite extensibility. Depending on the persistence length the polymer collapses 
into spherical (small persistence length) or toroidal (large persistence length) globule. It was shown that unfolding of 
both spherical and toroidal globule occurs with a decay in the reaction force in the tadpole regime (the tadpole's head 
is correspondingly spherical or toroidal) with the following drop when the globule completely unravels. In addition, 
Craig and Terentjev also considered the model of ordered globule which can describe unfolding of inhomogeneous 
globules. 

In our recent wortP, we have performed a self-consistent field (SCF) modeling of equilibrium unfolding of a ho- 
mopolymer globule subjected to extensional deformation, using Scheutjens-Fleer numerical SCF (SF-SCF) approaclf^^. 
This method has allowed us to obtain conformational and thermodynamic properties of the considered system (e.g. 
the force-extension curves) in a wide range of polymerization degrees and solvent qualities and to reveal a general 
picture of the globule deformation in the constant extension ensemble. It was found that the globule unfolding occurs 
in three stages (Figure [l]), in accordance with the prediction of Halperin and Zhulina''. Our simulations indicate that 
the transition between weakly stretched and tadpole states is continuos whereas the transition from the tadpole to 
the uniformly stretched chain occurs junpwise and is accompanied by the force drop, in accordane with prediction 
of Cooke and William^. On the other hand, we have shown that for short chains, under moderately poor solvent 
strength conditions, unfolding occurs continuously: upon an increase in the imposed end-to-end distance the extended 
globule retains a longitudinally uniform shape at any degree of deformation. 

The aim of the present paper is to develop a quantitative analytical theory of globule unfolding in the constant 
extension ensemble that explicitly takes into account the (poor) solvent quality and enables us to analyze in details 
effects of finite chain length. This theory will allow us to go beyond the capacity of the SCF modeling and to calculate 
the globule-tadpole and the tadpole-stretched chain transition points (i.e. to construct the phase diagram) and to 
find corresponding conformational characteristics, reaction forces and force jump at transition, in a wide range of 
polymerization degree N and solvent quality expressed via Flory-Huggins parameter x- In contrast to earlier theories 
our model explicitely accounts for the induced asymmetry of the globule in the phase coexistence regime (we remark 
that the prolate shape of the tadpole head was recently found by Grassberger and HsrP^ in their analysis of the Monte 
Carlo simulation data). Using the developed theory we perform a comprehensive analysis of the globule unfolding 
that includes (i) calculation of force-extension curves and their quantitative comparison with those obtained in SCF 
modeling"'^], (ii) a detailed study of the phase transition, construction of the system's phase diagram, calculation of 
globule characteristics in the transition points (iii) derivation of asymptotic analytical dependences for characteristics 
of deformed globule in different deformation regimes and in the transition points. 

It should be emphasized that both the theory we develop as well as the SCF approach used inU-are of the mean field 
type. It is well known that the mean-field approach neglects fluctuations around the ground state for the considered 
systeirP^. This approximation is justified for large systems except of the vicinity of the phase transition points, but 
may fail for "small" systems. 

Although the kinetic aspects of the globule unfolding lie beyond the scope of the present paper, let us note that in 
real unfolding experiments and in computer experiments the distance between chain ends is increased continuously 
at some constant finite velocity. The extension rate plays an important role and may determine the choice of the 
unfolding pathway (see, for exampl^l^. From this point of view, our analysis corresponds to an infinitely slow globule 
deformation such that at each end-to-end distance D the system has enough time to reach the equilibrium state. 
However, our theory may be insightful for the experiments in which the unfolding kinetics is studied because it allows 
to estimate the height of the barrier separating different states (different minima) on the free energy landscape. 

The rest of the paper is organized as follows. In the section "Model and general formalism" we introduce our model 
and derive free energies for weakly elongated globule, for strongly extended (unfolded) chain and for the microphase 
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segregated (tadpole) state. The results of the calculations, including force-extension curves, system's parameters 
in the transition points and phase diagrams are summarized in the sections "Numerical Results" and "Analytical 
Results" that is followed by Discussion and Conclusions. 



II. MODEL AND GENERAL FORMALISM 
A. Unperturbed and weakly extended globule 

1. Unperturbed globule 

Let us consider a flexible polymer chain comprising N monomer units, each of size a, immersed in a poor solvent. 
The solvent strength is charaterized by the Flory-Huggins parameter x ^ 0-5- 

Following Ref.^* we consider the polymer globule in the "volume approximation", that is the globule is envisioned 
as a "liquid droplet" of a constant number density of the monomer units ip 

where N is the number of monomers in the globule and V is its volume. The correlation length of the density 
fluctuations in the globule ^ ~ ip^^a^^ is assumed to be much smaller that the size of the globule Rgiobuie V^l^ 
that is the case either for large N or for sufficiently poor solvent strength conditions. The globule free energy can be 
represented in the form 

Fglohule =llN + jS, (2) 

where /i is the monomer chemical potential (free energy) in infinite globule (/i is negative) and 7 is the interfacial 
tension coefficient (7 is positive). Her e and below fc^T is taken as energetic unit. 

Moreover, as has been demonstratecP^E^, the conformational entropy of the chain in the globular state also scales 
proportionally to the globule surface area. Hence, 7 also comprises a contribution of the conformational entropy to 
the globule's free energy. 

Due to the interfacial tension, a free (unperturbed) globule has the shape of a sphere with the radius 



= (—) 



1/3 

(3) 



and interfacial area 5*0 — 4ttRq. 

Within volume approximation all the partial parameterd^ (p, /i, and 7 which characterize the globular state do not 
depend on the number of monomers units in the globule N and are solely determined by the solvent quality (i.e. by 
the Flory-Huggins parameter x), that is 

<p = vix) 
M = m(x) 

7 = -fix) 

The depen dence of ip and fi on x can be easily found in the framework of the Flory lattice model of polymer 
solutiorfi^E^. For a solution of chains with the polymerization degree N and the volume fraction v = (pa'^, the free 
energy per lattice site is given by 

F,,te = ^logu-H(l-w)log(l-u)-hxw(l-t;). (4) 

If we cancel the first term in Eq. Q corresponding to translational entropy of polymer chains in solution, we obtain 
the expression for the free energy of the solution comprising one infinitely long chain. 

The equilibrium volume fraction v of the monomer units in the globule (which is found in equilibrium with pure 
solvent) is determined from the condition of vanishing osmotic pressure tt 



^2 d_ fF.,, 
dv 



= 0. (5) 
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This leads to 

X = -i^^^--, (6) 

i.e. the desired dependence (p = ip{x) = o,^'^v{x) is obtained in "inverse" form (we remark that the translartional 
entropy is ehminated from the latter equation). 

On the other hand, /i amounts to the monomer free energy change when it is transferred from the pure solvent 
(dilute phase, which is taken as a reference state for calculation of the chemical potential) to the globular phase with 
the concentration ip. Therefore, first we re-write the free energy Q per monomer: 

Fmonomer = = log(l - w) + %( 1 - w) , (7) 



then 



or, usmg Eq. ^ 



1-v 

M = F„iono,ner(v) - F„,ononier{v ^ 0) = + H log(l - v) (8) 



/x = 2+— log(l-i;). (9) 

V 



Finding the dependence of 7 on % is a more difficult problem that has not a complete analytical solution yet. 
However, for moderately poor solvent (close to the coil-globule transition point) where the polymer-solvent interaction 
free energy can be represented in terms of virial expansion Fint.site — ^ ■ + ^ ■ v^, B and C are second and third 
virial coefficients, respectively, a closed form analytical solution for 7 can be obtained in the framework of the Lifshitz 
theory of polymer globule a^ 



Exact numerical prefactor in Eq. (10) was calculated bv Ushakova et. al^ (see Eq. (16) and the following discussion 
in Ref.""^^). In our case, expansion of the free energy Q at small v gives B = {1 — 2x) , C ~ 1 and, therefore. 

7a2 = ^(l-2x)^. (11) 

On the other hand, the values of 7 as well as those of f and fj, as functions of x can be found using the numerical 
SCF approach which was successfully applied to model the globule unfolding in our previous worlPi]. The method of 
finding the parameters ip, fi, and 7 are described in Appendix, the numerical values of p, /i, and 7 for x values ranging 
from X = 0-8 to 2 are given in Table 1 of Appendix. These values are also plotted in Figure [2] together with p{x), 
fj-{x) and 7(x) given by Eqs (|6|, ([9| and ( |Tl| ), respectively. We see a very good agreement of the numerical results 
for p and fj, with the Flory theory in the whole range of x, whereas for 7 an agreement is observed only at small x 
(i.e. where the virial expansion is applicable). Moreover, for strong solvent (x > 1-2) we see that 7 has a scaling on 



(1 — 2x) different than the quadratic one predicted by Eq. ( 11 ). Therefore, in the following numerical calculations we 
will use the values of p, fi and 7 given in Table 1. 



2. Weakly extended globule 

Let us assume that upon extension the globule undergoes only the shape deformation whereas its volume is con- 
served: The sphere transforms into a prolate uniaxial ellipsoid (often called spheroid), the major axis of the ellipsoid is 
equal to the given end-to-end distance D. We also assume that in accordance with the earlier SCF modeling resultJli' 
in a wide range of N and x such a deformation affects only the shape of the globule, the partial parameters {p, fj,, and 
7) remain unchanged. The monomer density in the deformed globule and, therefore, its volume are assumed to be 
constant and controlled by the local balance of attractive binary and repulsive ternary monomer-monomer interactions 
and do not change upon the globule deformation. The shape change leads to the increase of the globule surface area 
in Eq. ^ and, correspondingly, to the increase of its free energy. 

The volume and the surface area of the ellipsoid are 



(12) 
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and 



S 



(13) 



respectively, where h is the minor axis of the eUipsoid and a = yjl — \P- jD"^ is the eUipticity parameter. The asymmetry 
of the eUipsoid (major-to-minor axes ratio) is given by 



■ h \2Ro 



3/2 



(14) 



where we denoted by x the relative elongation of the globule with respect to its unperturbed size: x := D/{2Ro) > 1. 
Then the ellipticity is expressed as 



Inserting Eqs (13)-(15l into Eq. ^ we obtain the free energy of the prolate ellipsoidal globule: 

Fgiobuie = fJ.N + jS ^ nN + 7S'o g{x), 

where 

1 1 :e2 



s 



AiiRl 2x 2 



,3 _ 1 



(15) 



(16) 



(17) 



Monomer chemical potential in extended ellipsoidal globule fgiobuie £^nd the reaction force fgiobuie ^.re defined as 



^Fgiobuie 
fJ-globule = 7^ = M + 7 



and 



fgiobuie 



dF„io. 



— g(x) + So9{x) — 



dx 



-f^=lSog'{x)^^2nRojg'{x), 



(18) 



(19) 
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respectively. Hence, the restoring force is proportional to the derivative of g{x) which is equal to 

x{x^ - 4) 



1 



1 



3x 

c2 + 2(x3 - 1) ^ 2(a;3 - 1)3/^ 



arcsm 



(20) 



Figure [3] shows the plots oi g{x) and its derivative g'{x). 
At small extensions 

g{x)^l+lix~l)^ 
5 

g'ix)^^{x-l) 



X -1 « 1 



(21) 



i.e. the deformation bears elastic character, and the free energy change and the corresponding reaction force depend 
on the absolute value of the deformation 



Fglobule ^ IJ-N + -f 



So + ^{D-2Rof 



(22) 



47r 

fglobule — l{D — 2i?o) , 



^globule — ^ ~t~ 



3iV 



^(I?-2i?o)-^. 
15 ^ ^' N 



(23) 



(24) 



The function g'{x) passes through a maximum at a; 2.19416 and then decreases with extension. Asymptotic 
dependences at a; » 1 are g{x) ~ and g'{x) ~ The force decrease with extension indicates an absolute 

instability of the system in this range of extension. This behavior reflects the well known Rayleigh instabilitj^Sl of 
a liquid droplet under uniform uniaxial extension. A connectivity of the monomer units within a polymer chain 
suggests a microphase segregation within the extended globule. In addition to the compact globular phase, the second 
microphase should correspond to another deformation regime, namely, stretched conformation of the chain (D 3> 2i?o). 

All the formulae derived above are equally applied to the globular microphase in the two-phase regime. Such a 
"depleted" globule loses certain fraction of its monomer units which are transferred to the extended phase, and the 
number of monomer units in the globule is n < iV and its extension is equal to d < D. The values of n and d 
are determined by the phase equilibrium conditions. To use the expressions obtained in this section in the range of 
extension corresponding to microphase segregated globule one should simply change N to n and D to d. All partial 
characteristics (1^9, /i, and 7) depend neither on the number of segments in the globular phase, n nor on the deformation 
d. 
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B. Strongly stretched chain (unfolded globule) 



Now consider the limit of large D, D ^ 2i?o, where the globule is completely unfolded and the chain is strongly 
stretched, Figure [T]c. 

Let us assume that the probability of a contact between different monomer units (which are not nearest neighbors 
in the chain) is negligibly small. This permits to use for theoretical analysis the mode l of ideal chain. The first theory 
of stretching of the ideal polymer chain was elaborated more than 50 years ag cP^. The theory was developed for 
the constant force ensemble; correspondingly, the Gibbs free energy was calculated. A merit of this approach was the 
automatic account of finite chain extensibility which allows one to go beyond the range of applicability of Gaussian 
approximation for the chain elasticity. 

Let us take advantage of this approach and "switch" temporarily from the fixed stretching ensemble studied in 
the present paper to the fixed force ensemble. Using the principle of equivalence of thermodynamic ensembles we 
shall interpret the extension vs. force dependence D = D{f) as an inverse form of the force-extension dependence 
/ = f{D) and switch from the Gibbs free energy Gchain to the Helmholz free energy Fchain- 

For the following comparison with the SCF resultiP we model the chain as a random walk on a cylindrical lattice 
under the action of a force / directed along the z-axis. Let the probability to make a step either in r or in z direction 
is given by Ai, the probability of a step in "rz" -direction, i.e. simultaneously changing both r and z coordinates by 
±1 is A2 whereas the rest, Aq = 1 — 4Ai — 4A2 is the probability to change only the angular coordinate (f>. Then the 
partition function of a monomer unit is given by 

w = (Ao + 2Ai)eO + (Ai + 2A2)e^'^ + (Ai + 2A2)e-/'^ 

1 (25) 
= 1 - 2Ai - 4A2 + 2(Ai + 2A2) cosh (/a) = 1 + — [cosh {fa) ~ 1] 

2k 

where k :— l/(4Ai + 8A2) Then the partition function of the chain comprising N monomer units is 

Z=ll + ^[co8h{fa)-l]\ (26) 



2k 

The logarithm of the partition function gives us the Gibbs free energy, 

Gcha^n = - log Z = • log 1 1 + ^ [cOsh { fo) - 1] | (27) 

Once the partition function and the Gibbs free energy are known, the average chain extension D corresponding to 
the applied force / can be immediately found: 

„ dG chain ,^ sinh( fa) , ^ 

^^-^=^°- 2fc + cosh(/a)-l 

Eq. (28) can be considered as an implicit dependence of the force / on the extension D, Figure |4j 
The monomer chemical potential follows from the Gibbs free energy 

^■cha.n = G chain/ N = ^ | ^ + ^ ^'^^^^ ^ | (29) 

In order to return to the fixed extension ensemble, use the following relation between Gibbs and Helmholz free 
energies (i.e. perform the Legendre transform): 

Fchain — Gchain + D ■ f (30) 



relation between / and D given by Eq. ( 28 1 



The free energy and the chemical potential in Eqs ( 27 ) , ( 29 ) , and ( 30 1 can be expressed as a function of D using the 



A simple form of thermodynamic functions' dependences on the extension D can be obtained at /a <C 1, i.e. in the 



limit of weak deformations. Then from Eq. (28 1 we find 



from what follows 



D^'^ fa (31) 



fchain — 2fc • ^^2 ' (^2) 
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f^chain 



(33) 



Fch a 



(34) 



i.e the elastic free energy has a Gaussian form. 

Remarkably all the expression derived in this subsection contain a single parameter k which defines the conditions 
of a random lattice walk for the free chain. 

Note that with an increase in the degree of extension D/Na, the growth of both F and / calculated by exact 
formulae, Eqs (30) and (28), is more pronounced than that in the corresponding approximate expressions, Eqs (34) 
and (32 1. This is a manifestation of finite chain extensibility. At maximum extension, D/Na — > 1, the reaction force 
tends to infinity as /a ~ — log (l — ]^), see Figure [i] 

In the two-phase regime, the part of the chain consisting of {N ~ n) monomer units is stretched and coexists in 
equilibrium with the globular phase containing n monomer units. The expressions obtained in this subsection can be 
used for description of the stretched chain segment in the two-phase regime provided N is replaced by — n and D 
is replaced hy D — d. The values of n and d at given D, N, and x ^-re determined from the equilibrium conditions for 



the two-phase system, see Eq. (36) below 



C. Tadpole conformation 

In order to describe a microphase-segregated tadpole conformation with coexisting globular head and stretched tail. 
Figure [l] b, we assume that the globular head comprises n monomers and has a shape of prolate ellipsoid with the 
major axis length equal to d. Correspondingly, the tail consists of — n monomers and stretched at the distance 
D — d. Then the free energy of the tadpole conformation can be written as 

Ftadpole = Fgiobulein, d) + Fchain{N - n, D - d) (35) 

Note that we consider coexistence of only one globule with stretched tail(s) "attached" to the globule at one or 
at both sides, respectively. This minimizes the globule surface at given number of monomer units n in the globular 
head. 



The tadpole free energy Ffadpoie in Eq. (35) should be minimized with respect to n and d at each given value of D. 
It is straightforward to show that the minimization conditions are equivalent to equality of chemical potentials and 
reaction forces in globular and stretched phases, respectively 

lJ-globule{n, d) = ^J,chain{N - n, D - d) ^^^^ 
fglobulein, d) = fchatn{N -n,D - d). 
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FIG. 5: Force extension curves for the globule with A'^ = 200 at various values of x. Solid lines - theory, dashed lines - SCF 
modeling. 



Phase equilibrium conditions, Eq. (36 1 allow obtaining all the equilibrium characteristics of the system in the two- 
phase regime. Inserting the equilibrium values of n and d into Eq. ( 35 ) 
tadpole. 



gives the equilibrium free energy of the 
Comparison of the latter with the free energies of one-phase states, i.e. with Fgiobuie and Fchain, enables us 



to find the boundary values of extension Di and D2, which correspond to the transitions from the elongated globule 
to the tadpole and from the tadpole to the stretched chain conformation, respectively. The range of extensions 
Di < D < D2 corresponds to the thermodynamic stability of the tadpole conformation. 



III. NUMERICAL RESULTS 



In the general case, complete exact solution for the model introduced above cannot be obtained analytically, hence, 
in the present section, we perform a numerical analysis of our system and compare the results with those obtained in 
Ref.El using SF-SCF approach. 

The numerical values of the partial parameters ip, /i, and 7 are found from independent SCF modeling of free 
globules, see Appendix and Table 1. f and /i can be also calculated using the Flory theory, see Eqs ([6| and where 
the dependences are given in terms of polymer volume fraction v related to the concentration as v = ipa^. 

The numerical value of the parameter k that defines the law of lattice walk and, correspondingly, the elasticity of 
the open chain, is also chosen in the way that enables us direct comparison of theoretical and SCF modeling results. 
Namely, the nearest-neighbor and next-to-nearest neighbor step probabilities Ai and A2, respectively, were set 

iiP 

as follows: Ai = A2 = 1/9. This gives k — l/(4Ai -I- 8A2) — 3/4, and it is this value that is used in the numerical 
calculations whereas in analytical formulae k will be kept as an independent parameter. 



A. Force-extension curves and comparison with SCF results 



Force-extension curves calculated using Eqs (16)-(20), (27)-(30), (35), and (36) are shown in Figures[5]and[6]together 
with those obtained in ref.^^^ using SCF numerical approach. 

As it can be seen from Figures [5] and [6j the theoretical curves reproduce very well the shape of the force-extension 
curves obtained in SCF modeling. In both cases, force growth in the ranges of small and large D and a weakly 
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FIG. 6: Force extension curves for the globule with A*' = 500 at various values of x- Solid lines - theory, dashed lines - SCF 
modeling. 

decreasing quasi-plateau in the intermediate D range terminated by a force drop are observed. 

Moreover, the theory parameterization based on the SCF calculation of partial parameters of the free globule as a 
function of x. Appendix]^ leads to a very good quantitative agreement between the theory and the SCF modeling. 

Remarkably there are small divergences between the theoretical and the SCF modeling curves in the vicinity of 
the transition points between one- and two- phase state (globule-tadpole and tadpole - open chain, respectively). We 
anticipate that these divergences are related to the assumption of the constant density of the globular phase used in 
the theoretical model. We will return to this model applicability issue later. 

B. Phase diagram 

The good quantitative agreement between the force-extension curves obtained by means of the SCF modeling 
and the theoretical ones allows us to use it for studying other properties of the stretched globule that were only 
qualitatively studied irfSl. Figures [t] and [s] show typical deformation dependences of the number of monomer units, n, 
in the globular head and its asymmetry, S, respectively. One can see a progressive increase in the asymetry reflecting 
the extension of the globule in the one-phase regime. Then, in the two-phase regime, the number of monomer units, 
n, in the globular phase starts to decrease approximately linearly with D, at the same time the asymmetry S of the 
tadpoles head increases but its magnitude remains not too large. When the two-phase region comes abruptly to the 
end, the number of monomer units in the head n is still large. After that, as D further increases, the globular phase 
does not exist anymore. 

The obtained dependences of Di and D2 on N and x allow us to construct a phase diagram of the system, Figure |9] 
With a decrease in the chain length, the range of stability for the tadpole conformation narrows and Di{N), D2{N) 
curves for two transitions meet at a certain point {Ncr, Dcr)- This is a critical point for our system and in this point, 
i.e. for certain chain length N — Ncr, the range of extensions where the tadpole conformation corresponds to the 
global free energy minimum degenerates into a single point At this particular extension, D = Dcr, the ellipsoid, the 
tadpole, and the open chain conformation have the same free energy, thus coexisting in equilibrium. 

The value Ncr also implies the minimal chain length at which microphase segregation within the extended globule 
can occur. Ncr depends on the interaction parameter x- it is smaller for more solvophobic polymer. 

Other properties of the system in the transition points can also be found. For example. Figure [TO] shows the 
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FIG. 7: Number of monomer units in the globular phase in the extended globule with A'^ = 500 at various values of x- 



dependence of the number of monomer units, n, in the tadpole's head at the tadpole-open chain transition point on 
the overall polymerization degree N . We see that n is an increasing function of N and a decreasing function of x- 
The jump in the reaction force at tadpole-open chain transition is shown in Figure |11[ its magnitude decays upon 
an increase in N . The asymmetry of the globular head of the tadpole in the transition points, Figure [T2j decreases 
and tends to unity as N grows; in the second transition point the tadpole's head has obviously a more prolate shape 
compared to that in the first transition point. 



IV. ANALYTICAL RESULTS 



In the present section we derive closed analytical expressions for the force-extension dependence in all deformation 
regimes and for the globule's properties in the transition points. This derivation is based on approximate expressions 
for thermodynamic characteristics of the elongated (ellipsoidall) globule and those of the stretched chain given by Eqs 
([22])-([24| and Eqs ([32]) -([34]), respectively. 

In particular, we use quadratic approximation for the surface area of the prolate ellipsoida globule as a function 
of deformation. The applicability range of this approximation is restricted by the condition ^5^) ~ 1 = 2; — 1 = 

(52/3 — 1^1^ gee Eqs ( 17 ) , ( 20 1 , and ( 21 ) and Figure [3] As it can be seen from Figure [s] this condition are reasonably 
satisfied. 

The Gaussian approximation for the elastic free energy of the stretched chain, Eq. ( |22] ) is obtained from Eq. (27) 
in the limit of /a <C 1. Even though Figures [5] and |6] indicate that in the relevant range of deformations fa ^ 1 
we apply the approximations, (32 1- (33) in the following calculations keeping in mind that it underestimates the free 
energy and the force (Figure [4]) . 

Deriving expressions below, we will also assume that both the polymerization degree N and the number of monomer 
units in the globular phase n are large, i.e. iV > n 3> 1. Correspondingly, we will restrict ourselves in giving mostly 
the main in N or n term and next in TV or n correction term in the expressions below. 
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FIG. 9: Phase diagram of unfolded globule (limits of stability for the tadpole structure; first and second transition points Di 
and D2) at various values of x- Dashed lines show scaling according to Eqs (451 and (541. 
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FIG. 10: Number of monomers in the globular head of the tadpole in the unravelling transition point at various values of x. 
Dashed line shows scaling according to Eq. | [5T| ). 



A. Equilibrium characteristics of the deformed globule in the two-phase regime 



In the two-phase region the free energy of a globule in the tadpole conformation is easily obtained using Eqs (351, 

{D-df 



(22), and (34): 



Ftadpoie ^ fJ-n + Air^rl + (d - 2ro)^ + fc ■ 

5 [N — nja^ 



(37) 



where the first term corresponds to the volume contribution to the free energy of the tadpoles head, the second and 
the third terms describe its excess surface free energy, wherease the last term describes the free energy of the extended 
tail. 

Recall that n is the number of monomer units in the globular head of the tadpole, d is the head's main axis size 
and ro is the radius of the spherical globule containing n monomer units: 



ro 



/ Q \ 1/3 



(38) 



By minimizing the free energy ( 37 ) with respect to n and d (or, equivalently, applying the equilibrium conditions 
(36)), we obtain 



87r7rg Sttj o ^ '^o _ , (D-d) 



3n 



(d- 2ro) ■ — = -k- 
15 n ^ nya'^ 



471-7,, „ N „, D-d 
5 [N ~ n)a'^ 



(39) 



Let us find the long axis of the tadpole's head. By substituting {D — d)/{N — n) from the second equation of the 
system into the first one, one arrives at a quadratic equation for (d — 2ro). It gives the following solution for d: 
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FIG. 11: Jump in the reaction force (restoring force) in the unravelling transition point at various values of x- 
shows scaling according to Eq. ( 57 1 . 



Dashed line 



Then the tadpole's head asymmetry is given by 




(41) 



i.e. in the two-phase state the globular head is elongated, its longer axes is larger than the diameter of a spherical 
globule comprising the same number, n, of the monomer units. Remarkably, the second (correction) term in Eq. (40) 
is independent of n. According to the x-dependences of 7 and fi this term is a decreasing function of x- A relative 
elongation d/2ro and related to it asymmetry parameter 5 — {d/2ro)^^^ of the tadpole's globular head grow with an 
increase in D because of a decrease in the number of monomer units in the globular head n. Figure [Sj 
The n{D) dependence is obtained from the second of Eq. ([39| and reads 



5k{D - d) 
27r7(d - 2ro) 




d) 



(42) 



In the two-phase region, where D ^ d, the decrease of n as function of D is close to linear, the slope of this dependence 
is independent of N and decreases with an increase in x (see Table 1 in Appendix and Figure [?]) . 
The reaction force / is obtained from Eqs (391 and (40) as a function of n and is given by 



7 / 47r 



1/3 



(43) 



This equation shows that the reaction force decreases with extension: as D increases, the number of monomers in the 
globular phase, n, decreases and, therefore / decreases. Figures [5] [6| 

Note that the correction term in Eq. ( 43 1 growing with a decrease in n includes the interfacial tension coefficient 7 
in the combination with the other partial parameters of the globule: 'y\iJ.\~'^ip~^^^ 




FIG. 12: Head asymmetry in the transition points at various values of x- 



B. Transition points 



In order to find the values of deformation, D, corresponding to the boundaries (transition points) between different 
regimes, the free energies of possible states should be calculated and then compared; then the boundary values of D 
are found as intersection points of F{D) dependences. 

The transitions observed in the system are those between one of the one-phase states (prolate globule or extended 
chain) and the two-phase state (tadpole conformation). Therefore, the system of equations for determination of the 
transition point includes the two-phase state equilibrium conditions, Eq. (36), and the free energy equality condition 

^tadpole Fglohule -^tadpole 



Taking into account Eq. ( 35 1 the free energy equality condition reads 



Fglobule{n, d) + Fchain{N - U, D - d) = FgiobuleiN, D) 

Fgiobuie{n, d) + Fchai7i{N - n,D - d) = Fchain{N, D) 



(44) 



By solving together Eqs (36 1 and (44 1 one obtains the threshold value of the extension D as well as the vlues of n, 



d in the corresponding transition point. 



1. Globule - tadpole transition point 



The critical extension for the prolate globule-tadpole transition point is obtained quite easily. According to our 
model, at this point the new stretched phase appears. This phase has an infinitely small size, hence, we should take 
the limit n 



N and c? — >■ D in Eqs ( 35 1 and ( 36 ) This gives us the transition point 




7 / 47r 

R V3iv^ 



1/3N 



(45) 



and the corresponding reaction force in the transition point 
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The asymmetry of the elhpsoidal globule (see Eq. (14|) in the transition point is 

,3/2 /TTTj / n,„ \ 1/3 




(47) 



Eqs (45)-(47l describe dependences of /i, and at the transition point on partial characteristics of the globule, 
/i, and 7, and on the parameter of lattice walks k. According to Eq. ( |45| t he globule extension at the transition 
point, Di, scales in the main term as Z?i ^ i?o ~ N^^^ at large N (Figure M; the absolute value of the (negative) 
correction term decreases as a function of N. The restoring force / varies in the transition point continuously as a 
function of the leading term grows with x, the correction term scales as iV~^/^ and decays with N. The correction 



term in Eq. (46) includes the already mentioned combination or the partial parameters 'y\fj,\~ ip~ ' . The asymmetry 
of the globule in the transition point decreases as iV grows (Figure 12), the same is valid for the asymmetry of the 
globular head of the tadpole in the two-phase region. 



2. Tadpole - open chain transition point 



Now let us consider the second transition corresponding to unravelling of the globular head of the tadpole. This 
transition occurs, as it follows from the force-extension curves. Figures [5j [6j at substantial deformations. Assuming 
that D ^ d and equating the free energy of the tadpole (37) and that of the stratched ("open") chain (34) we obtain 
teh following equation: 



ini + Anjrl + ^^(d - 2rof + k ■ 



5 ' " N{N-n)a? 

The value of D can be expressed from the second equation of the system ( [SO] ) : 

27r7 







D 



5k 



{N -n)a^{d-2rQ). 



(48) 



(49) 



Taking into account that d — 2ro is given by Eq. ( 40 ) and neglecting the correction to the head surface energy (the 
3rd term in the Ihs of Eq. ([48|) we get 



/in + A-K-frl + {l - -^^ 



1 



7 / 47r 



1/3- 



Keeping only the dominating contributions of the highest order in n we obtain 

^3/4 



n : 



27\ 



3/4 / 4^ X 1/4 



3(^2 



(50) 



(51) 



i.e. the number of monomer units in the head of the tadpole at the transition point scales as iV^/^, see Figure 10 
similar to what has been found for a dry globule in xev^. The size of the tadpole's head, d, is easily found from 
Eqs (l40l) and (IsTl) 



27\ 



1/4 



27r7a 



from what the asymmetry of the head in the transition point follows 

5/4 / Q \ 3/4 



1 + 5y^k(p 



MV 

27 y 



iV-l/4. 



The extension D2 at the transition point is then obtained using Eqs (49), (51), and (52) and reads: 

3/4 / Q \ 1/4 



Do 



■ Na 



1 



3 1^27 



Ar~i/4 



(52) 



(53) 



(54) 
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The reaction force of the tadpole is obtained using Eqs (43) and (51 ) 

2^% 



ftai 



dpole 



1 



3/4 / 4^ X 1/4 



7V-1/4 



The reaction force of the stretched (" open" ) chain in the transition point is given by 

2kD 



hain 



2 VImI 



Hence, the force jump in the transition point is given by 

— ftadpole f chain 



3(^2 



and decays as N increases (Figure 10 1 



Ar-i/4 



(55) 



(56) 



(57) 



C. Critical point 

Using the formulas obtained above we can also find the critical point Ncr- It can be formally estimated using 
the expressions (45) and (54) for Di and Z?2 in the large globule limit (although the globule consisting of iV = Ncr 



monomers may hardly be considered as a "large" one). 



Nr. 



3/4 



(58) 



In the vicinity of the theta-point all partial parameters scale as certain power law functions of relative deviation from 



the theta-point, t = x/xe ~ 1- In particular, ip ^ ra ^, ^ t'^ 
Nt"^ ^ 1, that is the requirement of stability of the globular stati 



Hence, condition N 3> Ncr implies 



V. DISCUSSION 



A. Statistico-mechanical analogy 

Let us start with discussion of analogy between intra-molecular conformational transition and microphase coexis- 
tence in extended polymer globule and well-known classical problem of gas-to-liquid transition. 

The conformational rearrangements in the deformed globule are similar to the liquid-to-gas transition that occur 
with the van der Waals gas in the (V, T)-ensemble below critical temperature (T < Tcrit) upon increasing the system 
volume V . At small V , the system is in the liquid state which is the counterpart of the globular state. Figure [13] a. 
An increase in volume gives rise to the vapor phase coexisting in equilibrium with the liquid one, similar microphase 
coexistence is observed in the tadpole state of the globule. Figure 13 b. As V further increases, the molecules move 



from the liquid to the vapor phase and finally, when V is large, the system is completely in vapor (gas) state, in our 
system this role is played by the stretched chain state, Figure [T3|c. 

From the point of view of conformational changes, the analogy illustrated by Figure [13] is obvious. However, 
thermodynamic characteristics in our small (single polymer chain) system behave differently than those in macroscopic 
gas-liquid system. Comparison of p = p{y) and / — f{D) dependences. Figure 14 indicates a marked difference. 

In both systems, the pressure p and the reaction force / vary monotonically with an increase in V or D, respectively, 
in the one-phase states. (These changes however differ in sign because of different directions of reaction force and 
pressure). An essential difference between single polymer and macrosystem is observed in the phase coexistence 
regime. In the macrosystem [(V, T)- ensemble] the pressure p is constant at any system volume where two phases 
coexist, Figure[T4]b. In the extended globule the reaction force dependence on D is "anomalous": an increase in D in 
the phase coexistence regime is accompanied by a decrease in /, Figure [14] a. Although such a dependence is usually 
interpreted in statistical physics as a manifestation of instability, for the considered small (single polymer molecule) 
system this state cannot be avoided in the constant extension ensemble. In the range of D corresponding to phase 
coexistence the decrease in the force is monotonous but at the right edge it drops abruptly. The system therefore 
demonstrates a negative extensibility (a negative elastic modulus). 
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FIG. 13: Correspondence between states of the van der Waals gas in {V, T)-ensemble and conformations of stretched globule 
in D-ensemble. 
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FIG. 14: Globule force-extension curve obtained using SCF numerical approach for = 200, x ~ 1-4 in Ref.^^ (a) and 
schematic P — 1/ diagram for the van der Waals gas (b). 



B. Force-extension curves and phase diagram 



Let us discuss the obtained results on equilibrium globule stretching by comparing them with results of the SCF 
modeling reported in detail iiP. As it can be seen from Figures [s] and [6[ the quantitative theory is in a very good 
agreement with the SCF modeling results. 

Both approaches lead to the conclusion about existence of the three-state mechanism of extention of a homopolymer 
chain which is immersed into a poor solvent (% > 1/2) and has the globular conformation as an unperturbed state. 
Three states of the deformed globule correspond to a weakly elongated (ellipsoidal) globule and to a stretched chain at 
small and large respectively, and to coexistence of these two states within a single macromolecule at intermediate 
extensions. 

As follows from Figures [9] and Eq. (58 1, the system has a critical pointy i.e. there exists a minimal chain length 
Ncr below which, at < Ncr microphase segregation in the globule is not possible. Here N^r has the order of tens 
of monomer units &t x — 1-5 2.0 (the vloume fraction of the globule in this case is (pa^ ^ 1) and increases to 
approximately 200 monomers when x is decreased to 0.8 {i^a^ ^ 0.5). 

The three-stage mechanism is realized provided the condition N > Ncr is fulfilled. On the contrary, a small 
( "overcritical" ) globule upon extension is deformed "as a whole" , without intramolecular segregation. As follows from 
analytical expression (58), the value oi Ncr is independent of 7, i.e. it is determined only by the internal characteristics 
of the globule, but includes a dependence on fc - the characteristic of the free chain walk. 



Nc: 



3/4 



1/2 



Conformation of extended globule at these overcritical conditions was studied using SCF modelin^^^ [A^ = 200, x — 
0.8]. It was shown that in this case the assumption about a constant polymer density ip within the deformed globule 
cannot be fulfilled. Extension of a small overcritical globule causes not solely the shape deformation but also a volume 
change - correspondingly a decrease in cp and depletion (or "dissolution") of its core. Emergence of this state clearly 
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FIG. 15: Force vs degree of extension curves at various values of A'^ and x- Dashed lines show the position of the plateau at 
TV oo. 



shows its thermodynamic advantage in extended small globule at weakly poor solvent. At the same time the value 
of Ncr shows the limit of applicability of our model since the assumption about constant polymer density ("volume 
approximation") is no more valid at small N, where the width of the interface (globule's interfacial layer) becomes 
comparable to the radius of the globular core. This is also illustrated by the deviation of (p values found in the SCF 
modeling of free globules at smaller N (large 1/A^) from the best fit on Figure 18 in Appendix [A| This explains the 
effect observed in the SCF modeling. 

When the condition N > Ncr is fulfilled, the second phase corresponding to the stretched (open) chain emerges at 
relatively small extension of the spherical globule, D — 2i?o ~ N^, see Eq. (HSj). In the two-phase regime, the number 
of monomer units in the globular phase n decreases approximately linearly with an increase in D, Eq. (42). This 
leads to the growth of the positive term in ^giohuie, see Eq. p4| ) due to the surface effect, proportionally to 7. The 
reaction force / weakly decreases with an increase in _D, see Eq. ( |43[ ), Figures [5) and[6| With an increase in TV or/ and 
a deterioration of the solvent quality (i.e. increase in x), the range of stability of the two-phase state increases even 
in reduced coordinates D/N, Figure 15 the negative slope of the "quasi-plateau" diminishes, and the force-extension 
curve acquires the shape typical for systems undergoing a phase transition by passing via two-phase state (see Figure 
Tllb). 



For the chains of finite length the quasi-plateau in the force-extension curves, Figures [5] and [6] ends up with an 
abrupt drop of the reaction force. As it can be seen from analytical expressions, existence of the drop as well as the 
negative slope of the quasi-plateau is associated with the excess surface free energy of the globular phase. Positive 
contribution of the surface to the free energy of the two-phase (tadpole) state makes it at Z3 > D2 less advantageous 
than the one-phase state of the stretched open chain. According to Eq. (51 ), the minimal size of the globular head in 
the two-phase regime scales as n ^ N'^^'^ and as it follows from the results of the quantitative theory, this value is quite 
large. This value of n appreciably exceeds the number of monomers necessary to form a stable free globule. The 

latter value can be easily estimated from the condition that the free energy of the free globule, Fgiobuie = A*-^ + 7-^0' 
should be non-positive. This gives (with the account of Eq. (|3|): 



|pi|3(/j2- 



(59) 



Remarkably in the long chain limit {N ^1) the size of the minimal globule in the extended chain n does not depend 
on the elastic constant k of the open chain and its dependence on the solvent quality includes the same combination 
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of partial parameters as in the case of the free globule: 

3/4 



7 



iV3/4 , Nrmn ^ • (60) 



This combination of parameters also determines the quasi-plateau slope of the force-extension curve, Eq. (43). As it 
can be see from Table 1, as x increases, n and Nmin decrease. 

Note that on f{D) curves obtained in the SCF modeling (Figures 5, 6) some peculiarity is observed, namely, a 
little peak in the vicinity of Di where the new phase corresponding to the stretched chain emerges. The monomer 
chemical potential in the extended phase, ^chain, is independent on the total number of monomer units in this phase, 
but depends on the force /, see Eq. ([29])). Therefore, the results of the quantitative theory lead to continuity of the 
force at this transition point. An insufficient advantage of forming a tail (stretched chain phase) of very small size 
should lead to a certain widening of the stability range of the one-phase elongated globular state. It is also possible 
that the initial withdrawal of the stretched chain segment from the globule is accompanied by a change of the local 
shape of the globule surface. 

More visible is the difference between the results of the theory and the SCF modeling in the second transition 
point which occurs at strong deformation, i.e. in the vicinity of Z?2. According to the SCF results the transition 
from tadpole to uniformly stretched chain is shifted towards stronger deformations, as compared to prediction of the 
theory. 

This effect can be explained in the following way: At strong extensions (close to D2), the number of monomer units 
n in the globular head of the tadpole is small and, as a result, the deformation affects not only its shape, but also the 
density and other partial characteristics, similarly to that occcurs upon deformation of a small overcritical globule. 
As a consequence, force jump in the transition point observed in SCF caclulations turns out to be lower than that 
predcited by our theoretical model. 

C. Tadpole stability: about the role of the tail 

Let us recall that in the framework of the developed theory we study equilibrium characteristics of the systems. 
Therefore we can equally consider both globule unfolding by increasing D and globule refolding by decreasing D. 
In the latter case, the transition from the one-phase stretched chain state to the two-phase tadpole state occurs at 
extension D2 and corresponds to formation of a new globular phase. Because of the interfacial tension the globular 
phase formation in the extended chain becomes possible only by drawing a large number of monomer units n > rijnin 
into the globular phase. This jump-wise change in the number of monomer units in the stretched phase leads 
correspondingly to the jump in the reaction force at D2 on the force-extension curves. 

Analysis of the results of the theory and/or those of the SCF modeling shows another interesting effect. Consider 
the following mental experiment. Figure |16| Let us take a single polymer chain composed of A'^ = 500 monomer 
units and put it into a poor solvent with x = 2 where it collapses into a dense [fa^ ~ 0.96) spherical globule. Then 
we start an unfolding experiment: take the ends of the chain and increase the end-to-end distance D. According to 
the force-extension curve in Figure 6] and phase diagram in Figure [9] at D/a « 11.3 the mirophase segregated state 
appears, and the globule acquires the tadpole conformation. The tadpole can be extended up to I?/a « 297, Figure 
16 a, the globular head at this extension contains n « 112 monomers (more than 20 % of the total amount of polymer) 
and its major axis is equal to d/a « 7.3. This is the limit of the tadpole stability: a larger extension leads to the 
tadpole breakup and to the transition into the stretched state. Figure [T6| b. In other words, n w 112 is the smallest 
globular head that we can obtain by extending the globule for the given N and x- Consider the 200th monomer of 
the chain: since the tail is extended uniformly, the position D of this monomer can be found in a straightforward way: 
D/a = 7.3 -f 500-112 ■ (200—112) sa 73. Let us "pin" the 200th monomer at this position and cut (delete) the monomers 
from 201 to 500 inclusive. Figure [TB] a, c. Now we have a globule containing N — 200 monomers and extended at 



D/a = Ti in the tadpole conformation. Figure 16 c. According to our "cutting procedure" the equilibrium conditions, 
Eq. (36), will be fulfilled and no rearrangement in this "reduced tadpole" should occur. However, in contrast to the 
chain comprising 500 monomer units this tadpole can be further extended, from D/a > 73, up to D/a ~ 101. The 
latter elongation corresponds to the limit of stability of the globule with N = 200, where the head contains n 68 
monomers. Figure [T6| d. 

VI. CONCLUSIONS 

Motivated by recently obtained results of the SCF modeling^, we have developed an analytical theory of mechanical 
unfolding of a homopolymer globule in the constant extension ensemble, i.e. by imposing the distance D between the 
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FIG. 16: Scheme of mental experiment. 



ends of the macromolecule. Correspondingly, the (conjugated) observable in this case is the average reaction force /. 

Our approach consists in calculating and then comparing the free energies of the possible conformational states of 
the deformed globule. A weakly deformed globule is modeled by a prolate ellipsoid with constant density whereas 
at strong deformation it is represented as an "open" freely jointed chain. At moderate extensions, the co-existence 
of (weakly prolate) globular and (strongly) extended phases within a single macromolecule occur, i.e., the globule 
acquires a "tadpole" conformation with globular "head" and stretched "tail" . 

Intramolecular microphase coexistence in the tadpole conformation implies equality of the monomer unit chemical 
potentials in two microphases and force balance between the globular head and the stretched tail which inevitably 
leads to asymmetry of the head's shape 

The model includes several parameters, such as degree of polymerization N, solvent strength quantified by Flory- 
Huggins parameter Xi the mean polymer density in the unperturbed globule (p, interfacial tension coefficient 7, and 
monomer chemical potential in the globule fi. In fact, setting N and x determines the values of tp, 7, and fi. We 
have calculated the values of ip, 7, and n numerically in the "infinite globule approximation" by using SF-SCF 
numerical approach. The values of if and n can also be calculated using lattice Flory theory of polymer solutions; 
both approaches show very good agreem ent. The value of 7 can be found analytically for moderately poor solvent 
using Lifshitz theory for polymer globuleJ^^E^. 

Parameters of a random walk on the lattice for the freely jointed chain have been chosen similarly to those used 
in previous SCF numerical modelin^^. This allows us to compare force-extension curves obtained using SF-SCF 
modeling and in the framework of the developed theory. The results show not only a qualitative but also a very good 
quantitative agreement. This correspondence is better, the larger are x and/or N. 

Both the SCF calculation^^] and the quantitative analytical theory prove that upon an increase in the end-to- 
end distance D, three regimes of deformation succesively occur: at small deformation the globule acquires a prolate 
shape, the reaction force grows linearly with the deformation, at moderate deformations the globule is in the tadpole 
conformation, reaction force is weakly decreasing, then at certain extension the globular head unravels and the globule 
completely unfolds, the reaction force drops down and then grows again upon further extension. 

Furthermore, the analytical theory developed in the present paper makes it possible to go beyond the limits of the 
SCF calculations. First of all, it allows calculation of force-extension curves for large N, where the system size is 
large and numerical SCF calculations become very time and memory-consuming. Moreover, in the framework of the 
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developed theory it is easy to calculate prolate globule-tadpole and tadpole-open chain transition points and to find 
corresponding conformational characteristics, such as the asymmetry and the number of monomers in the globular 
head, reaction forces and force jump at transition, in a wide range of N and x- Assuming Gaussian elasticity for 
description of the stretched chain conformation we are able to obtain asymptotic analytical expressions for the force- 
extension curves in all the regimes of deformation as well as those for the globule characteristics in the transition 
points as a function of the chain length N and solvent strength x- 

Our analysis has shown that the system exhibits a critical point, i.e. there exists a minimal chain length N^r below 
which, at < Ncr, the intramolecular microphase segregation in the extended globule does not occur. The globule is 
deformed "as a whole", without intramolecular segregation but by progressive "dissolution" of its core . At the same 
time the value of Ncr indicates the limit of applicability of our model since the assumption about constant polymer 
density in the globule ("volume approximation") is no more valid at smaller N, when the width of the globule's 
intcrfacial layer becomes comparable with the radius of the globular core. 

Finally, we may conclude that the simple and physically transparent model developed in the present work adequately 
describes globule deformation in the constant extension ensemble. 



Appendix A: Appendix. Determination of ip, fi, and 7 by SCF approach 

Consider a free (unperturbed) globule. It has a spherical shape and its free energy is given by 

Fglobule — + jSgiobule (Al) 

If the polymer density within the globule is constant and equal to ip, then its radius and the surface area 

Rglobule = ( 4:5^ ) ' Sglobule = '^T^R^lobule = ^TT {-^^^j (A2) 



and 



or, per monomer unit. 



Fglobule = A*^' + 7 • 47r (A3) 



\47riy9 J 



^=M + 7-(^)''V-/^ (A4) 

Hence, by performing SCF calculations on the 2-gradient cylindrical lattic^^ for various N and plotting the depen- 
dence Fgiobuie/N VS. Figure 17 the values of ^ (point of intersection with F/N axis obtaining by extrapolation 
to = 0) and A ^ (367r/^2)i/3 

7 can be found. To determine the intcrfacial tension coefficient 7, polymer density 
within the globule, Lp should be known. 

It can be found with the aid of the axial distribution of monomer units n(z), (the number of monomer units per z 
axis unit): 

n{z)=Y,L{r)<j>{r,z). (A5) 



where L{r) is the number of lattice sites with coordinate (r, z) and 0(r, z) is the polymer volume fraction profile 
obtained in 2-gradient SCF calculations. 

In cylindrical coordinates, the surface of the spherical globule with the radius Rgiobuie is described by the equation 
Rglobule — r'^ + z"^ If the polymer density within the globule is constant, then n{z) = (p ■ irr'^a — ip ■ T^{R^iobuie ~ ^^)'^ = 
n-max — ■ TTz^a where 



is the peak value on the n{z) profile. If n„iax is known, then ip is easily found 



'P ■ T^RllobuleO- = • 71" -T— « (A6) 



1 fir?"^ 

^ = 9^AV^ 
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FIG. 17: Free energy of the unperturbed globule per one monomer as function of (symbols) and the corresponding 

linear fit (solid lines). 



TABLE I: Values of polymer concentration (p, monomer chemical potential ^, and interfacial tension coefficient 7 calculated 
for different x using SCF method 



X 


TV 






70^ 


0.8 


1000 - 


4- 3000 


0.54 


-0.10 


0.088 


1.0 


1000 - 


4- 3000 


0.70 


-0.23 


0.18 


1.2 


1000 - 


4- 2350 


0.80 


-0.37 


0.27 


1.4 


500 


' 1650 


0.87 


-0.54 


0.38 


1.6 


500 4 


- 1150 


0.92 


-0.71 


0.48 


1.8 


100 ~ 


- 1000 


0.94 


-0.89 


0.58 


2.0 


100 ^ 600 


0.96 


-1.08 


0.670 



The value of (p corresponding to the infinitely large globule can be found by extrapolating the dependence ip{N) vs 
1/N to — 0, Figure 18 The values of /i, 7, and (p obtained using the above described procedure for the values 
of X used in this work are shown Table [T] We find a very good agreement between the numerical results and those 
obtained in the framework of Flory theory, see Eqs ^ , ^ and Figure [2] 
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FIG. 18: Average polymer density in the unperturbed globule as function of 1/N and the corresponding extrapolation to 
iV oo. 
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